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MODULATION EQUATION FOR SPDES IN UNBOUNDED DOMAINS WITH 
SPACE-TIME WHITE NOISE - LINEAR THEORY 

LUIGI AMEDEO BIANCHI AND DIRK BLOMKER 


Abstract. We study the approximation of SPDEs on the whole real line near a change 
of stability via modulation or amplitude equations, which acts as a replacement for 
the lack of random invariant manifolds on extended domains. Due to the unbounded¬ 
ness of the underlying domain a whole band of infinitely many eigenfunctions changes 
stability. Thus we expect not only a slow motion in time, but also a slow spatial mod¬ 
ulation of the dominant modes, which is described by the modulation equation. 

As a first step towards a full theory of modulation equations for nonlinear SPDEs 
on unbounded domains, we focus, in the results presented here, on the linear theory 
for one particular example, the Swift-Hohenberg equation. These linear results are 
one of the key technical tools to carry over the deterministic approximation results to 
the stochastic case with additive forcing. One technical problem for establishing error 
estimates rises from the spatially translation invariant nature of space-time white noise 
on unbounded domains, which implies that at any time we can expect the error to be 
always very large somewhere in space. 


1. Introduction 

We study the approximation of stochastic partial differential equations (SPDEs) on 
unbounded domains near a change of stability of a trivial solution via modulation or 
amplitude equations. Due to the unboundedness of the underlying domain a whole 
infinite band (i.e., an interval) of eigenfunctions changes sign and therefore the trivial 
solution its stability. Thus neither the classical theory of invariant manifolds for PDEs 
nor the recently developed theory of random invariant manifolds lllSlISSlITH IHlTl can 
be applied. 

Modulation or amplitude equations are a replacement to overcome the lack of in¬ 
variant manifolds, and they serve as a universal normal form depending only on the 
type of bifurcation. Being widely used in the physics literature, they are a tool to 
describe the evolution of the amplitude of the dominating pattern changing stability, 
where close to bifurcation we expect not only a slow motion of the amplitude in time, 
but also a slow modulation in space due to the band of eigenvalues changing sign. 

For deterministic PDEs this theory is a well-established tool. See for example [lHl26l 
|36l I29II for classical references, and the detailed comments later in this section. But 
hardly anything is known for SPDEs on unbounded domains. 

As a starting point in this paper we consider the stochastic Swift-Hohenberg equa¬ 
tion [fT0l[22l. which is a reduced model for the first convective instability in the 
Rayleigh-Benard model and serves as one of the main examples in which pattern for¬ 
mation is studied. It is given as 

du 9 9 9 9 9/9 

( 1 ) — = -{l + d^fu + e^vu-u^ + e^l^E 

at 
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on the whole real line with space-time white noise As we want to allow for periodic 
patterns, we do not assume any decay condition of solutions at infinity. 

The theory of higher order parabolic stochastic partial differential equations (SPDEs) 
on unbounded domains with additive translation invariant noise like space-time white 
noise is not that well studied, while for the wave equation with multiplicative noise 
there are many recent publications (see for example 112511141 fTSlIlSin and even more 
recent ones for parabolic equations with very rough noise 1120112111 . 

In many cases parabolic equations with noise are studied subject to a spatial cut off 
or a decay condition at infinity. This is the case, for example, in 111711 . where the cut¬ 
off is both in the real space as well as in the Fourier space. Another example is II19II . 
In [13 and in a similar way in 11271 12811 . the authors consider -valued solutions, 
where an integral equation is consider instead of a PDF. The choice of trace class noise 
in these examples implies that we have an L^-valued Wiener processes and thus a 
decay condition at infinity, which in both cases leads to more regular solutions. 

If we were to consider decay at infinity, we conjecture we’d recover similar results 
but with a point-forcing in the amplitude equation, due to the rescaling in space, 
needed to obtain the modulation equation. 

The scaling of the equation involves small noise of order and small distance 

from bifurcation of order 0{e^v). Due to the closeness to bifurcation, we expect small 
solutions and slow dynamics in time. Moreover, a whole band of Fourier modes around 
wave-number fc± 1 changes sign close to /t = 0, and thus we expect the dynamics to be 
given by a slow modulation of the complex amplitude A of the dominant pattern : 

u(t,x) ^ eA[E^t, ex) • -I- c.c., 

where c.c. denotes the complex conjugate of the previous term. We expect such an 
estimate to hold on the slow time-scale with t = 0 [e~^). 

The noise is chosen in a way that in the limit e ^ 0 both noise and linear instability 
do influence the dynamics of the amplitude equation. If we scale differently, we would 
lose one of the effects. The choice of space-time white noise is mainly for simplicity 
of the analysis, in order to avoid further technical difficulties, as we expect space-time 
white noise to appear in the amplitude equation in many cases of coloured and thus 
smoother noise, thanks to the scaling limit. For a detailed discussion on coloured noise 
see where large but still bounded domains of order 0[1) were treated. Moreover, 
on bounded domains O, with fractional noise as in [|T| or a-stable noise, the scaling 
of the noise’s strength is different, but the result itself is similar. 

1.1. Previous results. In the equation without noise = 0) Mielke, Schneider & 
Kirrmann [123 showed (see also II36II or numerous other publications by the authors) 

(2) -I- vA - 3 [A[ ^A. 

For the stochastic equation Q on large but bounded domains of size 0[e~^) Blomker, 
Hairer & Pavliotis [|3 derived the stochastic amplitude equation 

(3) djA = 4d|A + vA-3\A\^A+7], 

on a bounded domain of order 0(1) with complex-valued space-time white noise, 
although ^ could have been quite regular in space. The idea of splitting the solution 
in a Gaussian and a more regular part, which we will use in our approximation result, 
was already present in this paper, but due to boundedness of the domain there were no 
problems with growth at infinity. See also Mielke, Schneider, & Ziegra II30II for large 
domains and no noise. 

Remark 1. On bounded domains the noise has to be of strength O(c^) to get inter¬ 
esting results. With that scaling Blomker, Maier-Paape, & Schneider [jdj showed that 
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the amplitude of the dominant mode is independent of space and derived a stochastic 
ordinary differential equation (SDE) in C given by 

djA = vA-3\A\^A+ p, 

where ^ is a complex-valued white noise in time. 

Remark 2. Spatially constant noise does not act directly on the dominant modes, and 
thus for noise of order the noise term would just disappear in the amplitude 

equation, and we’d only recover the deterministic one stated in ([^. 

If we increase the noise strength to be of order 0[e) and set the noise to be spatially 
independent as ^ = a[3, where p is the derivative of a Brownian motion in R, we 
obtain a time only white noise, and in that case additional terms in the amplitude 
equation arise due to nonlinear interaction of the noise with itself in Fourier space. 
Mohammed, Blbmker & Klepel II31II obtained in this case 

dj-A = 4d^A +vA+la^A-3\A\^A, 

which was already predicted by Hutt et.al. 11231124II by using a formal centre manifold 
reduction. 

1.2. Nonlinear vs Linear. In this paper we study the linear case, as the first step 
towards a full theory of modulation equations. Also in the already cited results on large 
domains by Blbmker, Hairer & Pavliotis [31 this is an essential step towards the full 
nonlinear result, which is somewhat separated from the remaining nonlinear estimates. 
Although here in weighted spaces the nonlinear estimate does not seem to be that 
straightforward as the nonlinearity is an unbounded operator. In [3, the authors used 
a splitting of the solution into a slightly more regular part in and a Gaussian part 
that was only bounded in C°, but allowed for better estimates due to its Gaussian 
nature. Here we focus on the Gaussian part only, but in contrast to [3 we face the 
additional problem that solutions and thus error terms are immediately unbounded in 
the spatial direction for |x| ^ oo. 

Let us first state the mild formulation of 0: 

(4) u(t) = e'^mo- f ds+ , 

Jo 

where the semigroup generated by = —(1 -I- 3^)^ -I- ve^ and the stochastic con¬ 
volution (t), which is the solution of the linear equation, are defined and discussed 
in more detail in later sections. 

For a result on modulation equations we need to compare this to the mild formula¬ 
tion of 0, given by 

(5) A(T] = - J dS -I- i^43^+riT), 

with the semigroup and the corresponding stochastic convolution 

defined in terms of a complex-valued Wiener process W. 

The main result is to show that 

u[t,x] — [A(^e^t, ex) ■ e'^ -I- c.c.] is small, 

and to do so, we have three key steps: 

• Nonlinearity - We need to show that we can control the difference between the 
nonlinear terms in 0 and 0. This should be similar although quite technical 
to the deterministic case and we do not treat this here. 
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• Initial Conditions - For the two terms containing the initial conditions Uq and 
Aq, we will see that they split in a more regular part that is treated by the 
known deterministic results and a less regular Gaussian part, which is not 
good enough to be treated by standard deterministic methods. We discuss 
these estimates in detail in the proof of Theorem |4.2[ 

• Stochastic Convolution - The difference between these terms is the main new 
result of this paper. Theorem |4.H It is the key estimate to prove a full approx¬ 
imation result for stochastic modulation equations on R. 

Let us remark that we focus on bounds in sufficiently good norms here. We might be 
able to give much simpler bounds in Lj^^-spaces, but then we would not be able to 
control in this norm the cubic —u^ later. Thus we focus on the supremum-norm, which 
is a good compromise. It does not require any order of spatial derivatives, but it still 
good enough to bound the nonlinearity. Unfortunately, in our weighted spaces, the 
nonlinearity is always an unbounded operator, so some more care will be needed here. 

1.3. Structure of the paper. In Section we introduce the stochastic convolution, 
and discuss its rescaling to the slow time-scale. We can already identify all error terms 
that need to be handled in the following sections. We do not follow the approach of 
Walsh II37II but the one of Da Prato & Zabzcyck llllll . using explicit series expansions 
for calculations. 

Before getting to the main results, in Sectionj^we present the key technical results 
for the stochastic convolution and the Gaussian initial conditions. We have an error 
estimate in spatially weighted C° spaces and its extension to estimates in space and 
time. The main assumptions are bounds on the Fourier kernel of the convolution 
operator, which are provided in the final three sections 

Section [^provides the main results of the paper. First we establish the approxima¬ 
tion result for the for the Ornstein-Uhlenbeck process (stochastic convolution), then we 
provide the key steps for the full attractivity and approximation results for the linear 
stochastic equation. 

Sections [^ - [^ provide, as already mentioned, the technical bounds on specific 
Fourier kernels that are necessary to apply the results of Section to the main re¬ 
sults. In Section [slwe present the bounds necessary in order to derive estimates in 
space. In Section]^ we provide the technical results necessary to obtain space-time 
estimates for the stochastic convolution. Finally, in Section we show the technical 
estimates that are necessary to treat the Gaussian part in the initial conditions. 

2. Stochastic convolution 

In this section we study the stochastic convolution introduced in First we intro¬ 
duce all spaces and definitions that we need for our analysis and then give the precise 
definition of the stochastic convolution. Finally we rescale it to the slow time-scale. 
This is necessary to identify the Wiener process driving the modulation equation, and 
provides the ansatz which error terms need to be bounded. 

2.1. Notation and Definition. In this part we present basic notation and definitions. 
We introduce all spaces used in the following and define a semigroup generated by our 
differential operator in terms of Fourier-multipliers. 

For some small y > 0 we define the norm 

||u||co = sup{(l + x2)-r/2Kx)|} 

and denote by C° the space of all locally continuous functions u :R —> R, such that 
||u||co < 00 . Analogously, we denote by C°j, the space of locally continuous functions 
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u:[0,T]x. 


such that the following norm is finite 


U\\cO = sup sup 

se[ 0 ,T]xeR 




Furthermore, we define the space by the norm 

lluik" = sup SUp{L“’'||u||cO([_i,i])}. 
^ se[0,r] ieN '■ ^ 

We also use the time independent version ° 


Lemma 2.1. The norms || • and || • \\cgo^ are equivalent for all j > 0 and all T > 0 
with constants depending only on y. 


Proof Let us remark that it is sufficient to verify the equivalence of || • H^o and 11-11.^0. 
First, for x e [—L,L], 

(1 + x2)-’'/2 ^ (1 + L2)-r/2 ^ ’i-fliL-r 


and thus 

sup (1 + x2)“’'/2|u(x)| ^ 

xe[-L,I] 

which easily implies the first bound. For the other bound just note that x e [—1,1] 
implies 

(1 +x2r»'^2|u(x)| ^ ^ ||li||.^o, 

and |x| e [L,L + 1] implies 

(l + X^r’^/^InCx)! ^ ^ 2’'||u||.^o. 

□ 


Note that for y < p we obviously have the following continuous embeddings 


^0 

"p,r 


^0 

-'r.T 


and 


C° c C°. 
p r 


Thus a bound on C° for any small y already provides bounds for all larger y. 

For small y > 0 these spaces are (up to a small e-dependent constant) almost in¬ 
variant under the rescaling x ex. The following result is straightforward to verify. 


Lemma 2.2. For any u e one has 

ceniulk;, ^ lluCeOlk", ^ C||u||.^o^ . 


The key observation is here that after the substitution z = xe one has 

e’'(l -F z^)-!'''^ < (1 + z2g-2)-r/2 < ^ ^ z2)-r/2 _ 


2.2. Local H^-spaces. We show that we can continuously embed (see below for 
the definition) into the weighted space C° for any positive p, as it is already embedded 
into the space C°(R) of continuous bounded functions which embeds into any C° with 

p > 0: 

||u||co = supj—|u(x)|} 

[equiv.] ^ sup ^ sup||u||c 0 ([_i,i]) = sup||iz||cO([i_i+i]) 

N3L>0 l€N L€Z 

[Sobolev] ^ C-sup||u||h 1 ([l,l+i]) = C-||u||hi [by definition]. 

L€Z 
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2.3. Semigroups and Green’s function. Here we recall well known facts about semi¬ 
groups in terms of Green’s functions and Fourier multipliers. Fix G^[x) to be the 
Green’s function (fundamental solution) associated to the differential operator 


if: 


-(1 + 


The semigroup generated by if is thus given as 




Gt(x-y)/(y)dy. 


We can write down G explicitly. Using the Fourier transform we immediately see that 
Gf = with gt(fc) = and thus 


( 6 ) 


Gt(x) 


g,ik)e''^^dk 

Jr 


1 


' 

^ikx 

Ju 


This is similar for the operator if^ = —(1 + + vs^. It has the fundamental 

solution: 

and via Fourier transform: 


G,t(x) 



+ ^ikx 


with kernel gg_t(fc) = e k^)^+ts^v^ 


2.4. Properties of the semigroup. The semigroup generated by if is a strongly 
continuous semigroup of linear operators II33II . Here we rely mainly on the explicit 
representation described in the section above. First we recall a bound on in L°°- 
topology and extend it to the weighted spaces. 

Using a result of []2 Lemma 2.1] we obtain that for all /3 > 0 there is a constant 
G > 0 such that 

ti/4|G,(x)|e(^+2t-'/bkl ^ c for all t e [0,1] 
and thus []9l Remark 2.2] 


(7) 


' 

|Gt(x)|e/’l’^ldx^G. 

JR 


With these bounds and some less optimal bounds for t ^ 1 from II31II . we immediately 
obtain the following Lemma: 


Lemma 2.3. For j e [0,1) there is a constant G > 0 such that for all t ^ 0 and all 


u^C° 


lle^-^'ullco ^ Gmax{l,t’'/2}||u||co . 

Proof First it is easy to see that we need to bound 

'' - l-F (x — z)^ A r/2 


I ^ _i_ ij^_ Z ) \ y/^ 

sup(l+x2)-r/2 |G,(x-y)Kl + yVdy = sup ( ) |G,(z)|dz 

I TO Ai€lR TO ^ 


For t e [0,1] the result now follows from Q and 

l + (x-z)2 


l-Fx^ 


^2(H-z^). 


For t > 1 we use that from II31II 


Gf(z)=t ^I^g,-ii 2 {t with sup sup |gT(?)K4-F ^ G . 

CeK Te(0,l) 
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Thus 


(l+z2)r/2|G^(z)|dz^C 


^ C 


^^^(4 + z^/t) Mz 

f (l + t?2)r/2 


4 +e 




□ 


Let us remark without proof that the restriction to 7 e (0,1) in the last step of the 
proof above does not seem to be necessary, as we could use any power of (4 + z^). 

2.5. Definition of cylindrical Wiener process. Let’s fix for the whole paper an ab¬ 
stract probability space (fi, on which all stochastic processes are defined. Fol¬ 
lowing mill we define: 

Definition 1. A standard cylindrical Wiener process W(t) is given by any orthonormal 
basis of L^(R,R) and any family of (real valued) i.i.d. standard Brownian 

motions such that 

( 8 ) W[t) = Y,Be[t)ei. 

teN 

Obviously, a cylindrical Wiener process is not an L^(R)-valued random variable. It 
is just defined in a larger space. For details see llllll . Moreover, it is characterized by 
being a Gaussian process such that for all u, v e L^(R,R) and all t,s ^ 0 

E{W(t),u) = 0 and E|{W(t),u)(W(s),v)| = min{t,s}{u,v) . 

A sometimes confusing fact is that for every fixed t the process {W(t,x)}j.g]g is a sta¬ 
tionary and thus translation invariant process, but the basis functions in which we 
expand are not at all translation invariant and might decay fast at infinity like Hermite 
functions. 

We also need the notion of a standard C-valued cylindrical process, which is not as 
standard: 

Definition 2. A complex-valued standard cylindrical Wiener process 'W{t) is given by 
any orthonormal basis of L^(R,C) and any family of C-valued i.i.d. standard 

Brownian motions such that 

(9) ir(t) = 2ft(t)e<. 

<eN 

Note that for a C-valued standard Brownian motion ft one has 
Eft(t)^ = 0 and E|ft(t)|^ = t. 

Example 1. If we take 2 independent copies j = 1,2 of real valued standard 
cylindrical Wiener processes in the sense of Definition then 

iV[t) := ^ -F iW®(t)] 

is a complex valued standard cylindrical Wiener process with same orthonormal basis 
of both L^(R,C) and L^(R,R) and with Brownian motion 

-F iB®] 

Pe = 


72 
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We will model space-time white noise always as the derivative of a standard cylin¬ 
drical Wiener process. There is also the equivalent approach of Walsh II37II by looking 
at the derivative of a Brownian sheet. But Dalang & Quer-Sardanyons II13II and also 
II18II showed that both formulations yield, up to taking the right versions of the pro¬ 
cesses, the same integrals. 


Remark 3 (Coloured Noise). As already mentioned previously, we only consider space- 
time white noise here in this paper. If we would consider a Q-Wiener process W, for 
example with Q given by the convolution against some function q, then we obtain for 
the corresponding noise ^ = d^W that it is a generalized centred Gaussian process with 
correlation 

E?(t,x)^(s,y) = 5(t -s)q(x - y). 

So we are still in the case of homogeneous noise which is translation invariant. Refer¬ 
ring to JB] the stochastic convolution is more regular in space, but we expect it to be 
still unbounded for every fixed t > 0. 

In []3 in the case of large, but bounded, domains the case of spatially smoother 
coloured noise is also treated, but the amplitude equation still displays space-time 
white noise. In order to avoid further technicalities in rescaling coloured noise, we do 
not address this issue in this paper. 


2.6. Stochastic Integrals. For a deterministic Hilbert-Schmidt-operator valued func¬ 
tion Jff e T^([S, T],IB(L^(R))), one can define the stochastic integral with respect to 
the real or complex cylindrical Wiener process W, which we can also expand in the 
basis Cf that were used to define the Wiener process in @) or 


rT 


^is)dW[s) = 


S 


dl3([s). 


By Ito-isometry we have 


E 


Jff(s)dW(s) 


rT 


E 

imJs 




ds = II J^ll 


2 


1.1. Stochastic convolution. For a cylindrical Wiener process W defined in @ and 
the semigroup generated by 5£ we can define the stochastic convolution 


( 10 ) 


w^(t,o 


' t 

Jo 



(eN-lo 


eO-O^eKOdftCs) 



Gt-5(--y)e«(y)dy dft(s). 


feN-'O JK 

This is by definition (see llllll ) the mild solution to the linear problem 

du = dt + dW, u(0) = 0. 


In this section we focus on W^(t,x), but the results are true for all other stochastic 
convolutions considered in this paper, like for example W.^g 2 [T,X) on the slow time- 
scale, with a complex Wiener process W. 

The following Lemma is well known in the setting of Walsh and easy to verify here. 
We give a brief sketch of a proof for completeness of presentation. 

Lemma 2.4. The stochastic convolution W^(t,x) in is for all t ^ 0 and x e R a 
weZZ defined real-valued Gaussian random variable with mean 0 and variance HgjU^ds. 
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Proof. It is easy to check that the series in ( [Tol ) is for all (t, x) e [0, T] x R a Cauchy- 
sequence in L^(f2,R) with 


E|W^(t,x)|2 


[Parceval] 


[Plancherel] 


E 

r f 


r f 


Gt-s[x - y)e([y)dy 
r f 


ds 


|Gt_,(x--)ll^ds = 


0 

r f 


iGJlMs 


IlgJI^ds < 00, 


where a straightforward calculation shows that the last term is finite. Gaussianity 
follows from the fact that Gaussians are closed under mean-square convergence. □ 


Remark 4. Lemma [2^ immediately implies that for all T > 0 and L > 0 

e LP(n X [0, T] X [-L, L]). 

With the methods of sectionjsjwe will see that e G°([0, T] x [—L, L])) with 

norm growing in L slower than any small power of L. 

The term behaves like the square root of logarithm at x ^ oo, a property we 
only state as a remark without proof. 


Remark 5. For all t > 0 

1 

lim - sup |W^(t,x)| e (0, oo) . 

ylogCL) xe[-L,L] 


This result has been obtained for Gaussian processes by Qualls and Watanabe Il34ll . 
The key idea required in proving this remark is just the fact that for fixed t > 0 we are 
considering a Gaussian field with constant variance (see Lemma 2.41. To conclude the 
proof one needs to compute and analyse the covariance function, which is a straight¬ 
forward but quite technical task. 


2.8. Rescaling. In this section we rescale the stochastic convolution to the slow 
time- and space-scale 

X = EX, T = e^t. 

We want to focus on the slow time scale and large spatial scale, forgetting about the 
local phenomena induced by fast oscillations and focusing only on the slow modulation 
of these fast pattern. 

Before we start, we rescale the Wiener process first, and show that such rescaled 

Wiener process is at least in law still the same cylindrical Wiener process. Note that = 
denotes equality in law. 

Lemma 2.5. For a real or complex valued standard cylindrical Wiener process W we have 
W[t,x) = 

teN 

where is an orthonormal basis in L^(R), {Pdeen ci family of (real or complex) 

Ltd. standard Brownian motions and we defined for some c e R; 

ef(X) = , /3f ^(T) = . 

Proof Check that is an orthonormal basis of L^(R). By the scaling properties of 
Brownian motion since ft(t) is a family of i.i.d. Brownian motions, a~^l3((a^t) is one, 
too. Thus by definition is a standard cylindrical Wiener process. □ 
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The full rescaling of the Wiener process is stated in the following Lemma. Here we 
do not start by but consider the final rescaled result, and assume that is indepen¬ 
dent of e, then we construct a proper c-dependent Wiener process for W^, such that 
the rescaling is true. 

Lemma 2.6 (Rescaling Lemma). Fix an orthonormal basis of L^(R,C) and 

{ftlteN n fatnily of C-valued Lid. standard Brownian motions. Then there exists an 
e-dependent standard real-valued cylindrical Wiener process such that 

e'^^^W^\Te~^,Xs~'^) = 

rT I- roo 

= — ^ e-iT-SX2+ksfe^ik[X-Y) dfce^(-y) dy dft -L c.c. 

Proof We define the complex Wiener process 

teN 

Define the projection onto the positive wave-numbers as the operator P"*" given by 

1 

P+/(x)=— dfc. 

2^ Jo 

We have the following identities, first for convolutions 

P+f ^g = jr-i [P+/g] = [/z[o,oo)g] =f* P^g, 

and then, using Plancherel, for scalar-products 

(P+u,v) = {u,P+v) = {P+u,P+v) . 

Let us now start with the rescaling 

1 r T p ^ 00 

2^ I Jo JkJ-i/e 

substituting fc' = (1 -I- ke) we get 

, _ rT r .. roo 

2^ I Jo Jk ^ Jo 

plugging in y = y/e and x = X/e yields 

1 r T ^ ^ 00 

= — 2 e-(^-®’"“'f^-'''’'e'''(^-^’dfce‘^e^(ye)dydft(S)-Lc.c. 

2^ ( Jo JrJo 

substituting now t = T/e^ and s = S/e^ gives 

^ r t r roo 

= — 2 e-('-'‘»-'^')'e'''(*-^)dfce'>e<(y£:)dydft(se2)-pc.c. 

2^ I Jo JrJo 

Using a rescaled complex Wiener process we have 

^ z't z' ^ OO 

= e— 2 e-P-Wi-kW^kix-y) dfce'^e^yc)dy dft(s) -L c.c. 

2^ ( Jo JrJo 


and, denoting the Green’s function by Gfx), this gives 



P+Gt_j(x - y)e'^ef[ye) dy dft(s) -I- c.c. 



MODULATION EQUATION FOR SPDE IN UNBOUNDED DOMAINS 


11 


We define now/^(y) = and again an orthonormal basis so 

that we obtain 


, 1/2 


s 

t 


r t 


* f([x)dP([s) + c.c. 


and by moving the P"*" it becomes 


=«‘"E 

t -1 


Gt_, * P+/tU) dft(s) + c.c. = ^ 

0 t 


r 1 




P feMdl3([s) + c.c. 


This is equal to E^^^W^\t,x') if we can verify that 

W'^'^\t,x) = ^P+//WftCO + c.c. 

( 

is a standard cylindrical Wiener process. Note that this would be trivial, in case 

i 

was a real valued Wiener process, as P+u + c.c. = u for any real-valued function. 

First we see that is a centred Gaussian process, we only need to check that the 
covariance-operator is the identity. 

E{w(®^(t),u){W^"^(t),v) 

= + c.c.,u)(P+fiPi[t) + C.C., v)] 

l 

= t2(F^,ii)(P+/t,v) + t2(P+/t,n)(^,v) 

t l 

= t2(P+/«,ii)(P+/t,v) + t2(P+/t,n)(P+/<,v) 

l l 

= t2{/^,P+u)(/^,P+v) + t2(/^,P+u)(/f,P+v) 

I I 


where we used that u and v are real in order to pull out the complex conjugate from 
the scalar-product. Now we use Parceval 


to obtain 

E(W^''^(t),u)(W^''^(t),v) 


l 

= t2(P+u,/<)(P+v,/^) + t2(P+u,/^)(P+v,/^) 

= t{P+u,P+v) + t{P+v,P+u) 

= t{P^u,v) + t{P+u,v) 

= t{P'''u,v) + t{P'^u,v) as V is real 

= t(u,v) 


as P'^'u + P^u = u for a real-valued function u. 


□ 


Now we observe that on the RHS in Lemma 1231 we almost have the stochastic 
convolution of an operator, if we suppose that ek is small and just had the innermost 
integral over R, instead that just (—1/c, -Foo). In that case we could have written 


(11) 


eW 4 g 2 (r,X)e'^ -F c.c. 
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What we want to do in the next sections is to give a bound on the error, the difference 
between what we have and what we want to use. Therefore, we provide first in Sec¬ 
tion!^ the technical Lemmas that reduce the task to calculations on Fourier-kernels of 
convolution operators. 


3. Holder Estimates eor general convolution integrals 

In ( [TT] ) all error terms are of the type where is a convolution 

operator written in terms of the Fourier-transformed kernel. For the initial conditions 
we also need bounds on H^A, with A being a Gaussian function. In this section we 
provide bounds on these objects in terms of norms of the kernel. 

First we need the following key Lemma for estimates: 


Lemma 3.1. Let W be a complex valued standard cylindrical Wiener process with or¬ 
thonormal basis in L^(]R,C) and a fatnily of C-valued Ltd. standard 

Brownian motions. 

Given a function /, its Green’s function H.^ = f [t, ■), and its corresponding con¬ 

volution operator = H.^* let us define 


HT) 


rT 

J^T_.,dW 

Jo 


expanded as 


rT 


Hr.AX-Y)ef(Y)dYd/3e(T). 


(eNvO Jm 

Then for all p > 1 and all y > ^ there is a constant C > 0 such that for all L ^ 1, 




where the L^(0, T,Wynorm of f is defined as 


2 


' r 

Jo Jr 


|/(S,fc)|2(|fc|2r + i)dfcdS. 


Let us remark that this is actually a slight abuse of notation, as we look at the 
H’^-norm of the kernel H^., which has Fourier transform /. 


Remark 6. Let us remark that similar estimates than the ones presented here were 
derived in 11141112II for the Green’s function of the stochastic wave equation and more 
regular noise. Moreover, we do not only need finiteness of the norms, but explicit 
bounds. Especially, the dependence of the constants on L. 


Remark 7. We will see that the conditions above on j are no problem, since thanks to 
Gaussianity we can do the estimates for p = 2 and then send p to infinity and hence 
we can choose y > 0 as small as we want. 


3.1 1 We proceed by steps. We start by using the fractional Sobolev 


Proof of Lemma 

embedding II35II and the explicit representation of the norm in 


= C 


11^(7, f-OII 




+ 


\HT,LX)-HT,LY)\p 


dXdy 


[- 1 , 1 ]" 


= CL 


-1 




[-L.Lp 


\HT,X)-HT,Y)\P 


dXdy 


with a fixed later such that y > a > 1/p. Let us remark that by Holder inequality it is 
sufficient to verify the claim only for sufficiently large p > 0. 
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For the next step we take the expectation and use the Gaussianity of 4>; 




( 12 ) 


+ CL“^+“P 


[-L,Lr 


Cp(,E\HT,X)-HT,Y)\^r'^ 
|X-y|i+«p 


dXdy], 


By means of this well-known trick we translated our problem from the generic p-th 
moment to the second moment only. 

We now proceed, bounding the second moments (first use Ito-isometry and substi¬ 
tute in time): 


teN-' 

rT 


(13) 


e\ht,x)-ht,y)\^ = Yj 

[ Farce val] 
[Plancherel] 


rT 


[H,iX-z)-H,iY-z)]e([z)dz 


dT 


\\H,[X--)-H,iY--)\\l.dr 


|/(T,fc)(e''^-e‘'^^)|2dfcdT 


0 

rT r 




0 J 


|/(T,fc)|2|fc|2r dkdT\X-Y\^L 


We used for the application of Plancherel, that for Hj. being the kernel of/(T, •) one 


has 


(14) 


H,(X-z)-H,(y-z) = 


/(T,fc)(e 


ik[X—z) _ ik(Y— 


"^)dfc 


/(T,fc)(e“-e'*^^)e-‘‘^" dfc 

& 

= ^[/(T,fc)(e''^-e‘''’^)] . 

Now we take ( [T^ and we plug it back into 


(15) ^ CL-iE||4.||^,^_,_,j 




\X-Y\TP 


We can compute 


|X - Y\TP 

[-L.Lh 




dXdY = C -LP^T-cd+i^ 


dXdY . 


so that ( [T^ becomes 

(16) E||$(r,-)|i;„ ^ CL-^Em[, + CLPr\\f\\[,^^^^. 

Now there is only one thing left: bound E||4>||^p. To do that, we take the same road 
followed above for the other term: 

•L fL 


( 17 ) 


E||4>(T,-)IIIpt_i,i]=E j \HT,X)\PdX 




(E\^[T,X)\^J/^ dX 


-L 


rL r 




( |/(T,fc)| 2 dfcdTj 


—L Jo Jl 


'dX^CLWfWl,^^^^, 
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and we can substitute that in obtaining the thesis for p large. To finish the proof 
for any p > 0, we use Holder inequality for q sufficiently large: 




□ 


Let us now extend to an estimate in space and time, where we first focus on a 
bounded domain in space and time. 


Lemma 3.2. Let $ be as in Lemma 3.1 and define for some y > 0 and T > 0 


rs 


II = sup S 
^ S€[o,r] 


- 2 r 


ll/(T,-)lli2 dT+ 


pR 


- 


+ sup (S - R) 

O^R^S^T 


||/(S-R+T,-)-/('t,-)lli2 dT . 


Then for all p > 1 such that Y > ^ there is a constant C > 0 such that for all L^ 1, 

« C ■ I"[ll/IIJ, + ll/lli/" 

Remark 8. Obviously, || • defines a norm, but we do not know whether it is equiva¬ 
lent to a known space. The first term has some similarities to Morrey-spaces, while the 
second one is an averaged Holder coefficient. But we do not need properties of that 
space, we just need to bound explicitly these norms. 

We recall the following Ito-isometry, which we already stated for the stochastic con¬ 
volution. 

Lemma 3.3. Let g be a square integrable function in space and time and W a complex 
standard cylindrical Wiener process. Then 


E 


^g(t)* dW(t) 


Ilg(t,-)lli2dt, 


where the stochastic integral is a function in x. 
Proof By direct calculation, 

h 

2 


E 


^g(r) * dW(r) 






r ^ 


£ J 
r ^ 


= Xi 1-^ ^g[r)*e(fdr 


i J a 


K -1 


g(r,x--), e<)^dr 


t 


= \W g(r,x--)lli 2 dr = ||g(r,-)lli 2 dr, 

J a J a 

where we used the Parceval theorem, translation invariance and Fourier isometry. □ 

Proof of Lemma [372] We use again fractional Sobolev spaces with y > a > 2/p and 
Gaussianity to obtain (with L^’-norm in space and time) 




rT rT rL 

[L-^\m[,+L-^+<^p 

Jo Jo J- 


^ CE 




-L J 

rT rT i*L 

Jo Jo J —L 




\Hs,x)-m,Y)\p 

Ks-Rf + 0(- y)2]i+“P/2 
Cp[E|$(S,X)-4>(R,y)|2]P/2 


dx dy ds dR 


KS-Rf + /X-Yf] 


21 l+ap/2 


dXdTdSdR 


]■ 
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Integrating equation (HZl) in time treats the L^’-norm in space and time, so we only 
need to look at the fourfold integral. Here we focus on the second moments in the 
integrand. Using ( [T^ , we obtain 

E|4>(S,X)-4>(R,y)|2 

^ 2E|4>(S,X) - 4>(S, Y)\^ + 2E|4>(S, Y) - 4>(R, Y)\^ 

^ 2\\f\\l,^^,^\X - y |2r + 2E|4.(S, Y) - m, Y)\^. 


Since the first term has already been treated before, we focus on the second one. If we 
can prove a bound by C||/||^ |S — R\^^, then we can easily finish the proof as in the 
previous Lemma 3.1 with sufficiently large p. 

Let’s assume without loss of generality S >R. 

/*S f*R 




Hs-t* dW[t) 


* dW(t) 




2 


^ 2E 

Hs_,*dW(t) 

JR 

+ 2E 



rR 


(Hs_,-HK_J*dW(t) 


So we have by Lemma [T3 


rs 


E|4>(s,y)-4>(R,y)|2 ^ 2 


rR 


= 2 


||/(S-t,-)||J.dt + 2 

0 

(*S—R (*R 

||/(t,-)||^.dt + 2 


||/(S-t,-)-/(R-t,-)ll?2cit 

||/(S-R+t,-)-/(t,-)lli2dt 


^ 2||/||2 -IS-Rpr. 


n 


The following Lemma states that the bounds previously obtained for fixed L > 0 
actually yield a bound in 0°^ j. 

Lemma 3.4. Fix T > 0 and y > 0 and let ube a random variable such that for all p > 0 
there is a constant C > 0 such that 

then for p > Vr 

L€N ^ 

Proof. Using Chebychev inequality yields 

P(||u||^o^^R) = p(3LgN : L-'niu||cO([o,T]x[-L,L]) ^ 

^ (ll^ilc“([o,r]x[-L,L]) ^ 

im 

Left 

□ 


3.1. Estimates for Gaussian initial conditions. In this section, we focus on technical 
results necessary to treat the term e^-^'Uo from the mild formulation, where the initial 
condition Uq[x) =A[ex)e'^ + c.c. is given by a modulated wave with centred Gaussian 
amplitude A. Due to Gaussianity, we need much less regularity of A than we would 
need in the deterministic case. 
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Lemma 3.5. Let A be a centred C-valued Gaussian with covariance operator Q given by 


a Fourier multiplier q(k) > 0. Let be as in Lemma 3.1 
Define 


= sup 
® se[o.r]J 


q(fc)|/(S,fc)P(|fc|2’' + l)dfc+ sup 

[ s,Re[o,r]J 


q(fc)-TT—-dfc. 


|S-R|2i 

Then for all p > 0, T > 0, and k > 0 there is a constant C > 0 such that 






Proof We know that A = for any orthonormal basis and standard 

C-valued Gaussians because any Gaussian can be written in terms of the co- 

variance operator Q and a cylindrical Gaussian Note that the symmetric oper¬ 

ator has Fourier multiplier q^^^. 

Using again fractional Sobolev spaces with y > a > 2/p and Gaussianity, it is suffi¬ 
cient to bound three second moments. First 


E| = E [Hsix - y) - Hsix - y)]A(y) dy 

' JR 

r 

= E 2 [HsiX - y) - HsiX - 

leN JS. 

= 2e|([Hs(X - --HsiX - 

leN 

= Y, I iHsiX - •) - Hs(X - •)], e£)L= 


;eN 
= llr)i/2 


= ||qi/2^-i[Hs(X--)-Hs(X--)]|| 


q(fc)|/(S,fc)P|fc|2rdfc-|X-y|2r, 


where we again used Parceval, Plancherel and in the final step ( [T^ . 
Secondly, we can proceed similar to the first case to obtain 


E| - M’rA[X)\'^ = 


q[k)\3F-\Hs-HiCl[k)\^dk 


^ sup q(fc)--—--dfc-|r-R|^y 

s.Re[o,r]jR IS-Rl^’' 


And finally, we can verify 


E\M’jAiiX)f = C 


qW\f[T,k)\^dk. 


□ 


4. Main Results 

This section provides the main results, where we focus on the new stochastic esti¬ 
mates and treat the deterministic estimates needed as an assumption. We give 

• Approximation Result for Stochastic Gonvolutions 

• Attractivity for deterministic initial conditions 

• Full linear Approximation Theorem 
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Let us first state the general formulation, and fix to be the proper rescaled 
Wiener process as in the rescaling Lemma 2.6 We consider a mild solution of the 
linear problem 

d,u = ^u + VE^u + e^^^d,w’^‘^\ u(0) = Uq , 

that is a function u which satisfies 

' t 

u(t) = e‘^UQ + VE^ + E^^^W^\t, x) . 

Jo 

In the whole paper, we need the following assumption: 

|v| ^ 1, i.e. V = ^(1). 

Defining ^g = ^ + vE^,we can rewrite the mild formulation in a slightly different 
fashion, that turns out to be easier to use: 

u(t) = + E^/^W^\tl 

We also have on the slow scale an approximating amplitude equation (AE) with a 
C-valued Wiener process and a solution A in the mild form: 

A(r) = + V f + W^giiT) 

Jo 

= + W^g2+,xn 

So now for the first step we assume initial conditions to be 0 and aim for the result for 
the stochastic convolutions only. 

4.1. Approximation Result for Stochastic Convolutions. Our goal in this section is 
||e^^^wjf^(t,x) - [W 4 g 2 _,_^(te^,xs)e'^ + c.c.Jllco is small 


for some small y > 0. Note that by Lemma 3.4 it is sufficient to provide bounds on 
moments of C°([0, T] x [—L,L])-norms. 

We 


Now we can throw in the definition of W^w, and rescale as in Lemma 


2.6 


obtain for some suitable convolution operator with corresponding kernel 

dW(s)e'^^'^ + c.c. = 


W^\Te-\Xe-^): 


^- 1/2 


271 


s 




(* 00 


g(r-s}[v-(2+te}^F] gifc(x-y) (-y ) + 


C.C. 


0 j] 

Thus 

gl/2^W(-yg-2^^g-l) _ W4g2+^[T,X)e^^^ 

rT' 

[^T-s - dw(s) 


3.1 


In view of Lemma 
kernel as follows 


0 


or Lemma 


^r!.ldW(S) 


3.2 


we 


define the convolution operator with 


y'^\x) = — 

^ 271 


r+oo 


^-z(2+kefk^+zv ikX 


r +00 


- 1 /e 


e"“ dfc- 

271 


g-T4fc2+TrgifcX 


1 r 

— e^'' 
271 L 


r +00 


[e 


-r(2+keyk^ _ -z4k^-\ ikX 


f-lle 


]e"^ dfc- 


,-z{2+ken\ikX 


]■ 
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Now we take the inverse Fourier Transform and we have (in view of Lemma [TT] ) 

/(T, fc) = e" 

Given 5 e (0,1) and using the symmetries of the integrals, it is sufficient to consider 


/, which is 

split into four pieces: 

(18) 

/(T,fc) = c-7(5/g.gg3(fc)e-^(2+;c.)^''^ 

(19) 


(20) 


(21) 



= fa+ fb+ fc+ fd- 


Here fa,fb> ^^d are the relatively simple terms that turn out to be small due to fast 
exponential decay. On the other hand, for we cannot take advantage of exponen¬ 
tially small terms, but we need to rely on the difference being small. 

We obtain the following main Theorem on the approximation of stochastic convo¬ 
lutions. 


Theorem 4.1. For all T > 0, for all k > 0, for all p > 0 and all sufficiently small y > 0 
there is a constant C > 0 such that 


P(||e^^^W^''^(t,x) - [W 4 g 2 _,_^,(te^,xe)e'^ + c.c.' 


^ ^ CeP 


for all E e (0, Eq). 

Proof Lemmajs^together with Lemmas 


5.1 


and 


-L, L]). Then the claim follows by using Lemma 3.4 


6.1 will provide bounds in C°([0, T] x 

□ 


Note that as we want to have the result for all small y > 0, we state here the result 
in C° instead of cP as in Lemma 

Y ^Y 

4.2. Attractivity for deterministic initial conditions. This section should motivate, 
why we assume in the full approximation result that the initial condition is a modulated 
wave that is split in a Gaussian with bounds on the covariance operator, that would 
not allow for the existence a derivative, and a more regular part in Hj^. 

We fix a time 

te = 

and try to approximate by a modulated wave. It turns out that we obtain a 

more regular part of the amplitude, and a Gaussian part, that is only in C°. 


3.4 


||e^«^«Uo + E^/^W^^\tJ - (Adet(fx)e- + A,,(ex)e- + c.c.)|| 




^ lle^'-^'Uo - (Adet(ex)e''' -F c.c.)||co -F J - (A,fEx)e^^ + c.c.)||co. 

For the first (deterministic) term we use the deterministic attractivity theorem, which 
we state just as an assumption. See [13^ for a result in which, together with the 
embedding ^ C° proved in Subsection 2.2 shows that the following assumption 
is at least true for all Uq 


Assumption 1. Suppose that for the initial condition Uq there is a smooth function 
such that 

— (Ajgt(£:x)e'^ -F c.c.)||(;o is small. 
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For the second, the stochastic term, we use our approximation result of Theorem 
l4.1l for the stochastic convolution to show that this is small. Thus we need to define 
^st as 


r ^ 

= W^ 32 ieh„X) ~ ^(o, ds) = ^(0,Q), 

Jo 

so we have a covariance operator Q with Fourier-symbol q given by 


q(fc): 


e ts 


ds : 


1 — e 




8fc2 


Direct calculation gives the following estimate on q(fc): there exists some Sg > 0 such 
that 


f 1 if |fc| ^ dn 1 

(22) q(fc)^c| ^ if |=Cmin{l, fc-} 

This will be needed as an assumption for the full approximation result. 


4.3. Full approximation. In this section, we will prove the general approximation 
result for initial conditions that are almost a modulated wave. To be more precise, we 
assume 

Assumption 2. Consider for some sufficiently small ^ > x > 0 the splitting 
Uo(^) = + C-C. + S^~'"E , 

where we assume for some sufficiently small k > y > 0 the following: 

\\E\\c» ^ C, ||AdetllH)„ ^ C, ~ ^( 0 ,Q) , 

where the Q is such that ( |22] ) holds. 

Let’s define for ease of notation Aq = Again, for ease of notation’s we 

define the solution u and the approximation as 

Uy^[t,x) = + W 4 g 2 (te^, ex)e'^ + c.c. 

For the deterministic approximation result we use the following Assumption. For a full 
deterministic approximation result in the space see for example [l36l . 

Assumption 3. Let’s define 

4et = sup [Ag,g,(ex)e-] - H^o. 

te[0,T„e-2] ^ 

We assume that this is a good approximation for the deterministic part, i.e. 

4^^ ^0 for s ^ 0 . 

Theorem 4.2 (Approximation). Under Assumptions^and^ for allK > 0 and y e (0, k) 
both sufficiently small and for all Tq> 0 and p > 1 there is a constant C > 0 such that 

pT sup ||u - u^llco ^ 4u + ^ 1 

[O.ToE-^] 

for all E e (0,1). 

We do not claim that the bound by e^^^ is optimal, but it seems that substantial 
improvements will be significantly more technical. 
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Proof. Using Assumption there are four terms in u — which we need to bound. 
There are three from the splitting in E, and and a fourth one from the differ¬ 
ence of the stochastic convolution. We proceed in several steps. 

First for the term with E we show by Lemma [23] for all t ^ 0 that 

^ Cs^-’^-r/^WEWco . 

r r 

Secondly, we approximate the difference of the two stochastic convolutions using The¬ 
orem |4A] Thus 

- [W^g2+,[T,X)e^^^ + c.c.]||co^ ^ 
with probability for all p > 1. 

Finally, for the remaining last two terms including Aq, we first use Assumption]^ for 
the term containing A^g^. 

For the final term containing A^j, we rewrite using a simple rescaling stated in 
Lemma 14.31 below 

e'^«(A,t(ex)e'3 = . 

Thus it remains to bound 


sup ||j^^A,t(ex)e'^||co < sup ||j^7.A,t(e-)llc» ^ sup ||j^7.A,t||co, 
re[o,ro] [o,r„] [o,r„] 


where we used Lemma 2.2 Now we proceed with Lemma 3.5 and Lemma 3.4 that 
gives 

E sup ||J#r4tllro < C||/||^ 


[o.r„] 


11.5?“ 


which is small as argued below in Section]^ Note that in view of Lemma 3.4 we need 
here the if “-norm for 7/2 and not 7, but we bound it for any 7 > 0 later anyway. 

□ 


For the terms involving the initial conditions we used the following rescaling lemma. 
Lemma 4.3. The following holds: 

^ .^j.A(ex)e‘^ 

with M’j having kernel 

f[T,i) = e- 

Proof We have that 


^ _ p-T{2+etfP+vT _ ^-4TP 


eT^-2^^[A[ex)e^^]= | -A((fc - l)/e)dfc 


e-TV\2+lsf-v]^(^^yisx 


Moreover, 


[g(43i-tv)T^](-gx).e^^ = 


d£ • e' 


□ 


Remark 9. Note that J#j. = M’j +Rt, where M’j has kernel /, the one we already 
introduced and studied in detail and an additional remainder Rj which has kernel 
f — f. In particular this kernel has only parts in the exponential tail, so the error is 
easily bounded. 

The error terms come from the fact that we do not cut in 0, but go further left to 
— 1/e. Also, we have contributions coming from the complex conjugate, but they do 
not cancel out. 
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5. Error estimates in space 


This section will provide the technical bounds on fa, ft, fr , an d defined in _ 

in the L^(H’')-norni. This is crucial for applying Lemma 3.2 in the proof of Theo¬ 
rem T 


4.1 These are all direct estimates that do not rely on any other result. 


Let us first remark that again we can bound separately the contribution with weight 
and with 1. Moreover, we can preform the computation for only, and then 
consider the case y = 0 to treat the 1. 

We use the following observations. For |v| ^ 1 and ek^ 5 we have 




and 


\fd\^Ce 


-T:4k^ 


Thus we obtain 


(23) 


l/a+/dl'|fc|'MfcdT^C 


0 


^ C 


^ c 


T /* 00 


,-8rf^2rdfcdT 


rT 

0 J 
rT 




-4Tk^ 


dk e 


5/e 


0 J 5/e 

-4z5^/e^ dT 

rTe- 


,-4T5^e-\-i2r+lV2 ^ ^ 


2-l-(2r+l)/2g-4<75^g2j^ 


(e O') 


0 


r 00 




-(2r+i)/2g-<745^ dcT = Ce^-^rs(.2r-i) ^ 


where, in order to be able to integrate in t, we must take — y — 1/2 > —1, that is 
y<l/2. 

By the same estimates we can bound the contribution of the term /j,. 


rT 


IA|2|fc|2rdfcdT = 


r-5/e 


rT 


lAl^lfcl^MfcdT^C 


-1/e 


„-8rf^2r 


5/e 


Now we turn to the complicated term f^. By using the mean value theorem we derive 

^- z { 2+kefe _ g-T4)c^ ^ [(2 -b kef - 4]fc2 

= —Te“'^^£:fcfc^(4+ fce) 

with E, taking value in [4fc^,4fc^(l + fcE/2)^], with the additional condition, given by 
the indicator function, that k e (—5/e, 5 /e). So the extremes of the interval for E, are 
actually 


ry 5>,2 52 52-, 

r 52 y 5>,2 52-, 




depending on k being negative or positive, respectively. We are interested in the abso¬ 
lute value of/, so we have, as t e [0, T], 


l/cK Z(-5/e.5A)e''^eTfc^(4-l-fce)e 
^ X{-5/e,5/e)^ 




^ X[-5/e,5/e)CeTk^e ^4*^' 


( 24 ) 
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Let’s now fix 0 < /x < 1/2 — 7 , in order to obtain 


/•T’ 


rT 


|/,|2|/c|2rdfcdT = 


/’S/e 


|/j2|/c|2rdfcdT 


-5/e 


rT 




/•5/e 


dfcdT 


-5/e 


(25) 


rT 




^ Ce^ 


^ Ce'" 


/•5/e 


0 

f * T /* 00 


0 

/•?■ 


fc2-'^fc4+2r+,.g-8T)c2 dfcdT 
^4+2r+/.g-8rf 




Remark 10. In the case v < 0 here and later many terms can be bounded indepen¬ 
dently of time, While for v > 0 our constants usually depends exponentially on T. 

So if we now put all contributions together we get the following bound: 

where we need y < 1/2 and /x -I- 2y < 1. 

Lemma 5.1. For all T > 0 and all jc > 0, there exist Eq > 0, pg > tmd C > 0 such that 




^ Ce 


l-K 


for all Y s (0, Yo) ^md all e e (0, eg)- 


In Lemma 


3.2 


6. Error estimates - 2: Time 

we provided a bound in terms of the norm ||/ 1 |^, which is defined as 


/•s 


II = sup S 

^ S€[0,T] 


-2r 


ll/(T,-)lli2 dT-l- 


/•« 


- Rl-2r 


-I- sup (S - R) 

O^R^S^T 


||/(S-R-l-T,-)-/(T,-)lli2 dT . 


Now we evaluate that explicitly. 

The first term is easily bounded as in Section We obtain the following bound for 
Y is close to 0: 


rs 


(26) 


sup S-^T 
se[o,r] 


ll/(T,-)ll?2dT ^ Ce^ _ 


The term did not appear in the estimates in Sectionjs] but is easily controlled. We 
comment on the proof in more detail below. When we are considering there was 
already a term in S showing, and we have from 


sup S ^ 1/2 / 2 / 2 ^,j .^^/2 gyp gl /2 /l /2 2 r^ 

se[o,r] Jo se[o,r] 

with /X < 1 to guarantee the integrability in 0 and /x = 1 — 4y to bound the supremum 
contribution by a constant and have the maximum /x possible, with y < 1 / 4 . For the 
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Other terms, we follow the estimates in ([23]), and we have 


sup S 

se[o,r] 


"*5 f* Ss 

ll/a,b,d(T,-)lli 2 dCT ^ sup S"^’'Ce 
0 se[o,r] jQ 


^-l/2g-45V 


dcj 




^ sup 

R€ [0,”|”Oo) 


a-i/2e-45Vdc7. 


Now R a dcr can be bounded with R for R > 1 and with R^^^ for 

0 ^ R ^ 1, and putting eveiything together we have ( |26| ). 

We can now move on to the second term. Using first that / is bounded by a constant 
and then the mean value theorem, we have that for some ^ between 0 and S — R and 
7] e (0,1), the second term is bounded by 


rR 


sup (S - R) 

O^R^S^T 


-2r 


||/(S - R + T, •) -/(t, •)lli2 dT 


^ sup [S-Ry-^-^C^ 


O^R^S^T 


[different;^] = sup (S— 


o<R<S:$r 


+ sup (S-R)^-"^C„ 


0^R:$S<T 


+ sup (S-R)^-"^C„ 


0^R:$S^T 


rR 

0 J: 

(*R I* 

0 Jm 

f^R 


0 Jr 


-^ifa+fd +fb +/c)lt=T+f 


dfcdx 


^(/a +/d)lt=T+C 


^ ^ C/b)lt=T+^ 


dfcdT (A+D) 

dfcdx (B) 


^ ^ C/c)! t=T+^ 


dfcdT . 


(C) 


We consider the three components separately. 

We start with B, and the same ideas will provide the bounds also for A and D. 


(B)= sup (S-R)’1-2’^C„ 


: sup (S-R)’1-2’'C^ 

O^R^S^T 


sup [S-Ry-^'^ 
o^R^s^r 

= c sup [s-Ry-^'^ 

OsiR^SsiT 

^ C sup (S - R)''“2>' 
o^R^s^r 


pR 


0 J 

f*R /*—5le 


^(/b)lt=T+c dfcdT 

.rU+C)v Ky _ (-2 + fce)2^2^g-(T+5X2+fcU";c^ |ri 


0 J-l/e 
rR ci/s 


(1 + dfcdT 


/•R 

0 d 
rR 


5/e 
r l/e 


(1 + dfce-’’4^^V'^^ dT 


5/e 


.^-(2T)+l)/2g-rj4T5Ve^ 


^ C sup (S - R)''-2>'(£:1-2’1 + e) 
o^R^s^r 


We need p ^ 2y, so we take p — 27 , and for the integrability in time we need j < 1 / 4 . 

The pieces (A) and (D), as anticipated, are bounded exactly in the same way, up to 
different constants, so we skip the details. 
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We are just left with the (C) part to estimate. Our best option to get rid of the 
difference in this case is to use the mean value theorem a second time, in k. We obtain 


rR 


(C) = C sup 

OsSR^S^T 


'5/e 

d 

J-5/e 

Q ^ C/c)! t=T+^ 


^ C sup (S - R)^ 

O^R^S^T 


— I gii(T+Ov 


5/s 


-5/e 


dfcdT 


y(^g-(T+0(2+fce)^ft^ _ g-(T+04fc^)^ 


(27) 


+ 


(^-(2 + fce)2fc2g-(T+?K2+fc£f ^ 4^2g-(T+C)4fc2 j 


dfcdT 


rR 


sup 

OsSR^S^T 


rS/s 


|^g-(T+5)(2+(te)2(c2 _ g-(T+C)4ft^^|?7_|_ 


-5/e 


+ 


- (2 + fc£:)2fc2e-('"+?)(2+(c^)2;c2 


dfcdT. 


As anticipated we have now two more instances of the Mean Value Theorem, for the 
functions e~“ and xe““ in the variable x, taking values in the interval 


I 




where the extrema of the interval might be switched due to the (additional) conditions 
on fc. But we can bound the size of the interval anyway (as already done previously in 
the space bounds): 


|f| 


4fc' 


(‘-(' 4 n 


= |fcpe|4+fceK5R|fc|^ 


We can now consider the two pieces of (C) (called (Cl) and (C2)) separately, one for 
each difference and write: 


rR 


(Cl)^C sup 


sup 

O^R^S^T 


r 




Ifcis^e-C^+^^P" dfcdT 


-5/e 


rs/e 


dfcdT 


where we took the value for p = 4fc^ (^1 — that would maximise the exponential, 

recalling that p e |^4fc^,4fc^ (l + • Moreover 0 < ^ < S —R. 

Now 


rR 


(Cl)^C sup (S-R)’’“2’'e''“''+'^ 
osSR^s^r 

rR 


r 5/e 


^2r,+/.g-rfr,4^4^ 


sup (S-R)’’"^’'e'" 

O^R^S^T 

^ C sup (S - ^ Ce'" 

O^R^S^T 


.j.-(M+l)/2-r, 


if —p/2 — p > — 1 / 2 , i.e. p < 1 — 2r], and we consider the optimal case, 17 = 2 y, so we 
can take p < 1 — 47 . This is slightly less then in the previous cases. 
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Now we do the same with the term (C2): 


(C2) ^ C sup (S - 


0 J 


-5/s 




sup 

O^R^S^T 


sup (S-R)''" 2 >'e'^ 

O^R^S^T 


(*R f* 5je 


0 J 


0 


/*R r 5/e 


0 


(1 + dfcdT 

^2„+Mg-r,rf dfcdT + 


■*R I* 5je 

+ C sup fc^^+'^e-^'^'^'dfcdT. 

o^R^s^r Jo Jo 

The first term is exactly what we had for (Cl), so we go on only with the second (C2.2) 


(C2.2)^C sup (S-R)''“2’'e'" 
o^R^s^r 


rR 


^-(5T)+(i+l)/2 


which is also very much alike (Cl), except from a slightly different exponent. In this 
case we need —2 < —5tj — /r — 1 ), i.e. /t < 1 — Srj, and by taking 17 = 2 y as before, we 
get /i < 1 — IO 7 , with Y small. So we have a final result analogous to the one for the 
first term ( |26l ), namely 


(28) 


sup (S-R)-^’' 

O^R^S^T 


rR 

11/(S - R + T, •) - /(T, Oil?. dT ^ crI-i'”'-''-. 


Putting together ( |26| and ( [28l ) we obtain the following bound. 

Lemma 6.1. For all T > 0 and for all k > 0, there exist Eq > 0, Yo> 0, and C > 0 such 
that 


Ik ^ 


for all Y s (0, Yo) tindfor all e e (0 , Rq)- 


7. Gaussian estimates 


In this section we provide the technical estimates that we need to apply Lemma 3.5 
in the proof of Theorem|4.2[ namely that the terms in 


||^„ = sup 
® se[o,r] 


q(fc)|/(S,fc)|^(|fc|^’' + l)dfc+ sup 

J S,Re[0,T]J 


q(fc) 


|/(S,fc)-/(R,fc)P 

|S-R|2t 


are small. We will consider the two suprema separately. 

Recall that we have a bound for q(fc) from ([22l): q(fc) ^ C min{l, k~^}. 


dk 


7.1. First supremum. We want to show for y e [0, VO that 


sup 

Se[0,T] 


q(fc)|/(S,fc)|V|fc| 2 r + l)dfc ^ C(r 2 + rVv . 


To prove it we use the form of /, that we know from equations ( [T^ , ( [T^ , ( [20l ) 
and As it was the case in Sections and we can bound the pieces f^, fj, 

and /j in the same way, and use a slightly different approach for f^. As in previous sec¬ 
tions, we can consider separately the terms with k^'^ and 1 , the second being a special 
case of the first one, when y = 0 . 
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In the first three cases we have: 


(A,B,D) ^ sup C 
se[o,r] 


q(fc)|e 


-4Sk^\2 




5/e 


00 

^ sup C 

Se[0,r] .5/e 


which is small, as long as 7 e [ 0 ,1/2). The remaining new terms from / — / are treated 
in a similar way. 

Finally it remains to treat f^, which requires some more care to treat, as the expo¬ 
nentials in the integrand cannot be bounded with a constant as it would result in a 
diverging integral. What we can do is to use the bound ( | 22 | ) on q(fc) and truncate in 
Sq to get rid of the singularity in 0. 

To be more precise using ( | 2 ^ we obtain 


rS/e 

q(fc)|fc|^’'l/c(S,fc)Pdfc^ q(fc)£:2s2|fc|6-t2rg-8S»c''^^ 

E J-5/e 

rSo rS/e 




k^+^'^dk + CE^S^ 




r 5/e 


^ Ce^S^ + Ce^S^ 


fc2+2r+(^dfc 


So 


^ Ce^S^ + 


provided that 0 ^ 27 -F /r < 1 . We can take p = 1/2 and we get for 7 e [ 0 ,1/4) 

•* 

q(fc)|fc|'"l/c(S,fc)|2dfc^C(e2 + eV.). 

JR 

Note that, in particular, the bound holds also when 7 = 0. 


7 . 2 . Second supremum. Also for the second supremum we consider two cases, as 
above, when splitting /. We want to prove that for 7 e [ 0 , V^) one has 


sup 

s,i?e[o,r] J 


q(fc) 


|/(S,/c)-/(R,/c)|- 

|S-R|2t 


dfc ^ -F • 


•Vb- 


We will use some of the estimates introduced in Section For the easy pieces, namely 
fa, fb and /d, we have 


r 00 


(A-FD,B)^ sup ClS-Rr^’' 

s,Re[o,r] 

^ sup C|S-R|^"^’' 
s,Re[o,r] 


q(fc)l/(S,fc)-/(R,fc)pdfc 


5/e 

f* 00 

5/e 

00 


/c-2|—/(t,fc)rd/c 
at 


^ sup ClS-Rl^-^’' fc2’J-2dfc 

S,Re[0,r] . 5/e 

^ sup C\S-R\'^-^^e^-^^, 
s,Re[o,r] 


which is small if 17 < 1/2 and if 17 ^ 27, so we can take q = 27. A similar estimate 
holds for the additional terms coming from f — f , since they are in the symmetric 
exponential tail. 
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There remains now just one case to check, what in Sectionj^was the term (C); we 
proceed in a way analogous to ( [27t . 

rS/s 


(C)^C sup \S-R\~^'^ 

S,R€[0,T] 


q(fc)|/(S,fc)-/(R,/c)|2d/c 


-5/e 
rS/e g 

sup IS-R^-^n q(fc)|—/(t,fc)rdfc 


s,Re[o,r] 


-5/e 

rS/e 


d t 


sup IS-Rl^-^^s'^ 
s,Re[o,r] 

r S/e 

[Assume 17 ^ 2y] ^ + k^^)dk 


q(fc)fc^^(l + fc2’’)dfc 


0 
r ^0 


^Ce^[ 


rS/e 


k^^[l + k^Vdk + 


dfc] 


[choose 17 = 1/5] ^ 

Remark 11. The choice of 17 = 1/5 might be optimised by choosing 5 -C 1. 
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